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Surface chemical processes on Si, particularly oxida- 
tion, are important for many technological applications. 
In the Si-based microelectronics industry, Si02 serves as 
gate dielectric. Therefore, the quality of the Si/Si02 
interface plays a significant role in the functionality of 
many Si based applications. Although most of the stud- 
ies of the growth of Si02 on Si have focused on ther- 
mal oxidation, native oxidation under ambient conditions 
is also important because during device production, H- 
passivated surfaces are often the starting point for se- 
quences of chemical and lithographic processes. There- 
fore, the durability of H-passivated surfaces is a critical 
factor in device processing. 

In addition to the technologically important (001) sur- 
face, oxidation of (111) oriented surfaces has also been 
studied. These studies have provided a vast amount of 
information regarding surface chemistry of this funda- 
mentally and technologically important chemical reac- 
tion. Availability of data from many different techniques 
makes (111) Si a convenient platform for studying sur- 
face chemical phenomena. For example IR absorption 
measures spectral changes in Si-H vibrations and relates 
them to chemical composition of the surface bonds. 
Ellipsometric measurements provide information about 
the thickness of the oxide, [l, Q X-ray photoemission 
spectroscopy |5|, l6| and photoelectron diffraction [7] pro- 
vide information about the suboxides at the Si/Si02 
interface. Reflectance difference spectroscopy measure- 
ments allow the number of oxidized Si layers to be 
counted during oxidation. [8] Each of these techniques is 
sensitive to different aspects of chemistry. 

A particularly powerful technique that has been used 
to study this interface is second harmonic generation 
(SHG). Since the dipole contribution is forbidden in 
centrosymmetric materials such as Si and amorphous 
materials such as Si02, the Si/Si02 interface provides 
the dominant contribution to the SHG signal in this 
system. Previously Daum et al. used spectrally re- 
solved SHG to study suboxides and their evolution dur- 
ing Si oxidation. [9, To] Downer and co-workers used 
SHG to probe excited hot carrier injection from Si-based 
films. [11| Tolk and co-workers showed that boron- induced 
interface charge traps can be studied using SHG.[T2l, [Tsf 
In these studies the SHG response is typically obtained as 



the sample is rotated, or the sample is aligned so that the 
plane of incidence is always parallel to a certain crystal 
orientation. Similar measurements performed in ambient 
conditions on (111) Si have provided information about 
early oxidation kinetics and passivity of H-terminated 
surfaces. p^fTBj 

We recently employed SHG anisotropy (SHGA) mea- 
surements to study bond-specific oxidation of Si. [l6|, [I^ 
By analyzing these data with the bond-charge model of 
nonlinear optics, we retrieve chemical kinetics character- 
istic of the different bonds. [18] We discovered that ten- 
sile and compressive strain speed up and slow down, re- 
spectively, the oxidation kinetics of different bond direc- 
tions, allowing control of in-plane surface chemistry dur- 
ing oxidation. [19| In addition, we found that the excess 
electrons in a heavily n-type sample change oxidation 
kinetics of the different surface bonds dramatically. [20| 
Here, we extend these studies to (111) surfaces of p-type 
Si wafers. We also performed atomic force microscopy 
(AFM) measurements, which provide additional infor- 
mation via changes in surface morphology. We find that 
p-type doping makes these surfaces very reactive to oxi- 
dation. The up-bonds are particularly affected, oxidizing 
very quickly and catalyzing the formation of the growing 
oxide layer. Despite the differences in chemical kinetics, 
the oxide terminates at the same thickness seen on n-type 
material. 

We investigated three (111) Si surfaces on boron-doped 
wafers with hole concentrations of 1.3 x 10^^, 3.6 x 
10^^, and 5.1 x 10^^ cm~^. These values were veri- 
fied by four-point-probe experiments. The samples were 
Hyrdogen-passivated by the procedure described in Ref. 
16 and then installed in an ambient atomic force micro- 
scope within 3 minutes of passivation and imaged in non- 
contact mode as a function of time. 

AFM experiments indicate significant surface reactiv- 
ity in air. Figure 1 shows selections from the time depen- 
dent sequence of AFM images for lightly (upper panels) 
and heavily (lower panels) p-doped samples. In these se- 
quences, the bright protrusions indicate the effect of sur- 
face oxidation as described previously. [21] The effect of a 
high density of holes is evident by comparison. In panel 
al and bl, the initial H-passivated surfaces exhibit flat 
terraces separated by monatomic steps. For the lightly 
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doped sample, oxidation takes longer than three hours. 
On the other hand, the heavily doped material oxidizes 
relatively rapidly. Most of the surface is covered with a 
structure consistent with a typical natural oxide in less 
than one hour. A high concentration of holes obviously 
impacts surface reactivity significantly. 

To gain further insight, we applied SHG to study the 
effect of p-type doping on the oxidation of the different 
surface bonds. The excitation beam was generated by a 
Ti-sapphire oscillator. It consisted of p-polarized 100 fs 
pulses centered at 806 nm arriving at a repetition rate of 
70 MHz. The p-polarized SHG response was detected by 
a photomultiplier at every 1° as the sample was rotated 
by 360°. 

Figures 2a and 2b show SHGA data for samples with 
hole concentrations of 3.6 x 10^^ and 5.1 x 10^^ cm~^ 
respectively. Both are normalized to the level where the 
SHGA responses no longer change with time. According 
to the bond model[18], the three major features at az- 
imuths 120°, 0°, and -120° arise when the alignment of 
the exciting field is roughly parallel to one of three back 
bonds. As reported previously [16], oxidation increases 
the magnitude of the SHGA signal at these azimuths by 
increasing bond asymmetry. A comparison of the two 
figures illustrates a significant dependence on carrier con- 
centration. The three features of the heavily doped sam- 
ple (Fig. 2b) rise almost linearly to their maximum value 
then remain constant thereafter. In contrast, the same 
three features of the moderately doped sample (Fig. 2a) 
reaches a higher maximum value and then decreases. We 
also note the existence of a 20 min incubation period for 
oxidation of the moderately doped sample. 

To obtain a more quantitative understanding of the re- 
action dynamics of oxidation, we investigated the time 
evolution of the hyperpolarizabilities for up and back 
bonds in more detail. The hyperpolarizabilities are deter- 
mined by least-squares fitting Eq. (26) of ref. 18 to the 
SHGA data, as described in more detail there. Figure 
3 summarizes the results for all three samples, thereby 
highlighting the differences that occur for different car- 
rier concentrations. With increasing concentration the 
curves shift to earlier times. At the highest hole concen- 
tration investigated here, both the incubation period and 
the peak near 2h have vanished completely. This shows 
that the up bond is now oxidizing as rapidly as the back 
bonds. 

We can make these observations quantitative by least- 
squares fitting simple exponential functions to the data. 
An expansion of the first 40 min of the data of Fig. 2a 
is shown in Fig. 4. The green curves show the expo- 
nential fits. The incubation period, if present, is fitted 
separately. The time constants for the incubation and 
rise periods are summarized in Table I. 

This analysis shows that the initial oxidation of p-type 



Si is dominated by the reaction dynamics of the up 
bonds and that holes facilitate the oxidation. For carrier 
concentrations of p = ~ 10^^ cm~^ as well as 10^^ cm~^ 
we see that the incubation periods for the back bonds 

TABLE L Parameters yielding the best fits of single exponen- 
tial functions to the data of Fig. 4. 





time constants [min] 




up bonds 


back bonds 


p [cm-3] 


incubation 


rise 


incubation 


rise 


1.3 X 10^^ 


-5.2±0.5 


33=b2 


-23=b2 


15=b2 


3.6 X 10^^ 


-1.2±0.3 


28±2 


-22±2 


13±2 


5.1 X 10^^ 


fast 


-22±4 


nonexistent 


-54±4 



are identical. However, the reactivity of the up bonds 
clearly changes, highlighting the dominant infiuence of 
holes during the first step of oxidation, where the H of 
the Si-H up bond is replaced by -OH. Here, the up bonds 
get through their incubation phase faster, the smaller 
negative incubation time constants of the up bonds com- 
pared to the back bonds indicate that up bonds reach 
saturation quicker than the back bonds. But their rise 
time constants are considerably slower than that of the 
back bonds once the back bonds start to oxidize as time 
evolves and -OH caps become available in the second 
step of oxidation. At this time the back bonds become 
more reactive and dominate the oxidation process with 
twice the rate of the up bonds. This is refiected in the 
rise times for these bonds. Positive time constants refer 
to new processes becoming active, the smaller these 
numbers get the faster the process is becoming active. 
For back bonds we obtain small positive time constants 
showing the dominance of these bonds once they become 
active. The extreme case is the p = 5.1 x 10^^ cm~^ 
sample. Here the incubation appears to be completely 
missing and the evolution is a simple exponential rise 
with a negative time constant, suggesting that the 
oxidation process is rising to a saturation level. Our 
ellipsometric data (not shown) show that for both n- and 
p-type samples oxidation continues until the thickness 
reaches a certain value, that is independent of carrier 
type and concentration |20|. These results are consistent 
with corrosion studies of Si and Ge in aqueous solutions, 
where the availability of holes which weakens the surface 
bonds is the rate- limiting step.[22l, |23| 

In conclusion, our AFM and SHGA data show that 
presence of holes in p-type material significantly in- 
creases the oxidation rates of H-terminated (111) Si 
samples, with the up bonds oxidizing at least as fast as 
the back bonds, and probably faster. 

The authors acknowledge D. E. A spnes for fruitful dis- 
cussions and J. Reynolds for the resistivity measurements 
of our samples. 
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I. FIGURE CAPTIONS 

FIG. 1. AFM scans of boron doped (111) Si samples 
for Oh, Ih, and 3h of oxidation in air (al)-(a3) carrier 
concentration p=1.3 x 10^^ cm~^ (bl)-(b3) carrier con- 
centration p=5.1 X 10^^ cm~^ 

FIG. 2. (a) Evolution of SHGA during air exposure 
of an H-terminated (111) Si sample with p=3.6 x 10^^ 
cm~^ (b) as (a) but for p=5.1 x 10^^ cm~^ 

FIG. 3. (a) Evolution of the average hyperpolarizabil- 
ities of the up bonds for boron-doped (111) Si samples 
with p = 1.3 X 10l^ 3.6 x 10^^ and 5.1 x 10^^ cm-^ . 
(b) As (a), but for the Si back bonds 

FIG. 4. (a) first 40 min of Fig 3a; (b) first 40 min of 
Fig 3b. The green curves show the exponential fits. 
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